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ABSTRACT 


Nnttleil  Mtkodi  for  the  analysis  of  the  elastic-plastic  fracture 
problem  os lag  a  special  finite  element  technique  are  presented.  A  brief 
description  of  some  ooacepts  in  elastic-plastic  fracture  mechanics  and  of  the 
finite  elwest  method  is  followed  by  the  formulation  procedure  of  the 
stiffness  matrix  using  eight-noded  quadrilateral  isoparametric  elements. 

After  a  terse  discussion  of  the  initial  stress  method,  the  procedure  of 
computation  is  extended  in  the  analysis  by  using  an  incremental  load 
proeess.  The  sixe  and  the  shape  of  the  plastic  zone  of  a  center  crack 
specimen  is  investigated. 

The  problem  of  eraek  growth  is  also  studied  using  the  moving  crack  tip 
element  procedure.  Results  are  presented  in  graphical  form. 

It  was  found  that  the  moving  finite  element  technique  can  be  modified  to 
solve  the  crack  branching  problems. 
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Introduction 

Fracture  mechanics  has  in  recent  years  become  an  independent  discipline 
that  deals  with  the  determination  of  the  conditions  under  which  machine  or 
strnctnral  elements  attain  nncontr oil able  failure  by  crack  propagation.  This 
topic  has  been  a  challenging  subject  in  the  area  of  solid  mechanics;  its 
historical  development  can  be  traced  back  to  the  early  part  of  the  century. 
A  knowledge  of  this  subject  can  assist  the  engineer  or  designer  to  safeguard 
machine  elements  against  catastrophic  failure. 

Two  fracture  criterions  of  failure  which  have  played  the  major  and 
fundamental  roles  in  the  study  of  fracture  mechanics  are  based  on  the 
concepts  of  energy  balance  and  stress  intensity  factors.  In  1921  and  1924, 
Griffith  [1,2]  introduced  a  fracture  criterion  of  failure  which  would  take 
into  account  the  effect  of  minute  flaws  and  small  cracks  in  material.  His 
theory  basically  involves  the  energy  balance,  which  states  that  the  energy 
consumed  in  creating  a  new  fracture  surface  within  solid  body  must  be 
supplied  from  the  release  of  the  elastic  strain  energy  and  from  the  work  done 
by  the  external  load.  Although  this  theory  has  since  been  proven  inadequate, 
it  still  provided  the  basis  for  newer  theories  developed  in  the  field  of 
fracture  mechanics. 

An  alternative  approach  was  later  proposed  by  Irwin  [3,4].  Rather  than 
follow  Griffith's  theory,  he  concentrated  on  a  crack  tip  region  which  is 
small  in  comparison  with  the  body  as  a  whole  but  sufficiently  large  with 
respect  to  atomic  dimensions  for  him  to  be  reasonably  comfortable  with  the 
application  of  linear  elasticity  theories.  He  stated  that  fracture 
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initiation  ocenrs  whan  the  intensity  of  stress  in  the  near  crack  tip  region 
reaches  a  critical  valne.  He  subsequently  introduced  the  three  nodes  of 
crack  extension  and  the  corresponding  stress  intensity  factors  Kjj,  and 
*III#  v*1*6*1  Provide  a  measure  of  the  amplitude  of  the  stress  field  in  the 
near  crack  tip  region. 

The  study  of  the  elastic-plastic  fracture  problems  are  the  main  concern 
of  this  report.  The  discussion  which  follows  will  be  restricted  to  Node  I 
(in-plane  tensile  mode)  crack  propagation.  The  problem  from  a  designer's 
viewpoint  is  to  prevent  brittle  fracture.  For  this  purpose,  the  state  of 
stress  and  strain  close  to  the  crack  tip  should  be  accurately  known.  The 
common  approach,  known  as  Linear  Elastic  Fracture  Mechanics  (LEFM).  is  to 
assume  linearly  elastic  material  behavior  around  the  crack  tip.  Since  it  is 
known  that  there  is  always  some  amount  of  plastio  deformation  at  the  crack 
tip.  the  LEFM  approach  can  provide  satisfactory  answers  only  for  highly 
brittle  materials  and  under  such  loading  condition:  which  cause  a  negligibly 
small  plastic  zone  at  the  crack  tip  as  compared  to  the  crack  length  and  the 
thickness  of  the  specimen.  Since  a  considerable  number  of  engineering 
materials  are  relatively  ductile,  and  since  the  thiekness  of  the  specimen 
cannot  always  be  considered  large,  the  plastic  deformation  at  the  crack  tip 
cannot  always,  be  ignored.  The  significance  of  the  size  of  the  plastic  zone 
ahead  of  the  crack  tip  as  one  of  the  parsmeters  which  govern  the  growth  of  a 
srr.:!:  ras  first  emphasized  by  Irwin  [3,4,9].  Later  researchers,  such  as  Rice 
[5,8,10],  Hutchinson  [6],  Cherepanov  [7],  and  others,  also  incorporated  the 
effect  of  the  shape  and  size  of  plastic  zone  in  their  proposed  fracture 
criteria. 

One  of  the  challenging  problems  is  to  investigate  the  <hape  and  the  size 
of  the  plastie  zone,  whioh  would  depend  on  a. ...  la!  properties,  specimen 


geometry,  loading  and  boundary  conditions.  A  considerable  amount  of  research 
in  this  direction  had  been  reported  in  the  literature.  A  good  bibliography 
on  the  subject  of  elasto-plastic  stress  analysis  can  be  found  in  references 
[101  and  [361. 

In  recent  developments  of  fracture  mechanics,  the  use  of  the  finite 
element  method  has  been  quite  extensive  in  both  elastic  and  elasto-plastic 
analyses.  A  number  of  special  finite  element  techniques  have  been  developed 
[16,17,18],  by  using  the  displacement  finite  element  method.  These  special 
elements  contain  a  singularity.  The  disadvantages  of  these  elements  are  that 
the  solutions  do  not  always  converge  and  that  they  take  too  much  computing 
time,  especially  for  elastic-plastic  fracture  problems. 

From  the  above  considerations  and  practical  viewpoint,  Barsoum  [11] 
formulated  the  eight-noded  quadrilateral  isoparametric  element.  The 
singularity  in  the  element  can  be  achieved  by  placing  the  mid-sided  node  near 
the  crack  tip  at  the  quarter  point.  It  has  been  shown  that  such  elements  in 
their  non-singular  formulation  satisfy  the  essential  convergence  criteria, 
namely,  continuity  of  displacements,  compatibility,  constant  strain  modes  and 
rigid  body  motion  modes.  Therefore,  their  convergence  characteristics  and 
the  high  accuracy  of  the  results  in  a  reasonably-sized  finite  element  mesh 
make  their  application  in  the  modeling  of  the  elastic-plastic  fracture 


problem  very  docile. 


SECTION  II 


Finite  Element  Analysis 


8  Concents 

fundamental  concept  of  the  finite  element  method  is  to  construct  a 
model  composed  of  a  set  of  piece-wise  continuous  functions  defined 
Lnite  number  of  closed  sub-regions.  The  sub-regions  are  called 
ts"  which  are  connected  to  each  other  at  their  common  nodes,  and 
rely  approximate  the  shape  of  the  domain.  The  generalized 
sent  components  of  these  nodes  are  the  basic  unknowns  in  the 
i. 

the  column  vector  {q}  be  the  generalized  displacement  components  at 
ss;  and  column  vector  IXJ}  be  the  state  of  displacement  for  the 
The  two  parameters  can  be  related  and  expressed  in  matrix  form  as: 

-  [NlCqJ  (1) 

I  is  the  shape  function,  which  is  a  function  of  special  coordinates 
to  be  evaluated  at  each  of  the  nodal  points.  The  notations  {  }  and 
:e  a  column  vector  and  a  matrix,  respectively. 

>ting  the  generalized  strain  vector  as  (e),  the  strain  components  can 
expressed  in  terms  of  the  displacements  through  the  compatibility 
is: 


-  IBJCD} 


(2) 
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where  [SI  is  s  partial  differential  operator  matrix  which  can  be  written  as: 


Through  the  constitutive  relations  of  stress  and  strain,  the  generalized 
stress  vector  {0}  can  be  written,  in  terms  of  the  strain  vector  {e},  as: 

{0}  =  [D] (e)  (4) 

where  [D]  is  the  material  matrix  which  is  given  below 
a)  for  plane  stress 


b)  for  plane  strain 
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where  E, v  are  the  modulus  of  elasticity  and  Poisson's  ratio,  respectively. 


2.  Quadrilateral  Isoparametric  Element 


The  formulation  of  the  stiffness  matrix  of  the  eight-noded  isoparametric 
element  is  well  documented  [28].  Through  the  following  transformation,  the 
stiffness  matrix  is  found  as  follows: 


y 


Figure  1.  Quadrilateral  Isoparametric  Element 


8 

x  -  Z  N.  (R,S)x. 


where  N.  ere  the  shepe  functions  corresponding  to  the  node  i,  whose 
coordinates  are  (x^,  y.).  In  the  transformed  R-S  coordinates  system, 

■=+lf or  the  corner  nodes  and  zero  for  the  mid-side  nodes. 


R2S2 

Ni(R,S)  -  [(l+RRi)(l+SSi)  -  (l-R2)(l+SSi)  -  (1-S2)(1+RR  )] 

S2  r2 

+  (l-R^d+SS^l-R2)  +  (1-S2)  (l+RRi)  (1-S2)-^ 


It  is  quite  easy  to  understand  that  if  the  boundaries  of  each  element 

2  2 

vary  from  -1  to  +1  in  the  R-S  coordinate  system,  the  terms  (1-R  )  and  (1-S  ) 
cannot  coexist  on  any  boundary.  The  interpolation  function  for  an  eight-node 
isoparametric  element  can  thus  be  written  in  the  following  form: 


-  -  ^(l+R) (1+S) (1-R-S) 

n2(r,s) 

-  -  ^(1-R)  (1+S) (1+R-S) 

n3(r,s) 

-  -  ^(1-R)  (1-S) (1+R+S) 

N4(R,S) 

-  -  i(l+R)£-S)  (1-R+S) 

N5(R,S) 

-  |(1+S)  (1-R2) 

N6(R,S) 

-  |(  1-R)  (1-S2) 

n7(r,s) 

-  ^(1-S)(1-R2) 

Ng(R,S) 

-  |(1+R)(1-S2) 
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Hie  displacements  are  interpolated  by: 

8 

u  -  Zni(R,S)u1 
i-1 
8 

v  -  Zhi(R,S)vi 
1*1 

The  strain-displacement  relations  can  be  written  as: 
/“i> 

<e>  -  (BJ  {  j 


where  [B]  is 


i 


(10) 


(ID 


(12) 


(13) 


where  [J]  is  the  Jacobian  matrix  and  is  in  the  following  form: 


[J]  - 


3x 

*L 

3R 

3R 

3x 

iz 

3S 

3S 

(14) 
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where  [D]  is  the  materiel  matrix  which  is  a  function  of  mechanical  properties 
of  the  material.  The  element  stiffness  matrix  CK]  is  given  as 


[K] 


(4 
•  -1 


[B]T[D] [Bjdet i J |  dR  dS 


(19) 


Since  the  form  of  the  shape  function  Nj(R,S)  in  all  elements  are  polynomials, 

3N.  3N 

the  derivatives,  —gg,  are  non-singular.  From  equations  (11),  (12),  and 

(13).  the  strain  can  be  written  in  the  following  form: 


{e} 


[J]_1[B'(R,S)] 


(20) 


Therefore,  from  the  above  equation,  the  singularity  can  be  obtained  by 
requiring  that  the  Jacobian  matrix  [J]  be  singular  at  the  crack  tip.  In 
other  words,  the  determinant  of  the  Jacobian  det  j jj  vanishes  at  the  crack 
tip.  This  can  be  achieved  by  placing  the  aid-side  mode  at  the  quarter  points 
of  the  sides. 


3.  fluilii  gftiai  Arwad  ttf  Jrtyfc  lifl 

In  order  to  obtain  the  singular  element  matrices  to  be  used  around  the 
crack  tip,  [s]  and  (a]  must  be  singular.  The  singularity  can  be  achieved  by 
placing  the  mid-side  node  of  any  side  at  the  quarter  point.  Figure  2a  shows 
the  two-dimensional,  eight-noded  quadrilateral  quarter  point  element  and 
Figure  2b,  the  six-noded  triangular  isoparametric  elements  with  the  mid-side 
nodes  near  the  crack  tip  at  the  quarter  point. 
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One  of  the  most  important  features  of  these  elements  is  that  they 
satisfy  the  necessary  requirements  for  convergence  [11]  in  their  singular 
form  as  well  as  in  their  non-singular  form.  They  possess  rigid  body  motion, 
constant  strain  modes,  interelement  compatibility,  and  continuity  of 
displacements,  in  contrast  with  many  other  special  crack  tip  elements. 

2  5  1 


Figure  2a.  Eight-Noded  Quadrilateral  Quarter  Point  Element 


Figure  2b.  Six-Noded  Triangular  Quarter  Point  Element 


It  is  probably  true  ia  saying  that  the  elastic-plastic  behavior  of  the 
material  played  a  dominating  role  in  recent  developments  of  fracture 
mechanics.  One  of  the  major  objectives  of  this  report  is  to  provide  a 
simplified  treatment  of  a  very  complex  problem  so  that  with  a  minimum  of 
programing  effort,  a  reasonable  solution  can  be  obtained. 

The  practical  interest  in  solving  the  elastic-plastic  fractnre  problem 
may  be  in  the  prediction  of  displacements  stresses  and  strains  at  various 
stages  of  loading,  the  size  of  the  plastic  zone,  and  the  growth  of  crack. 

The  method,  which  serves  the  purpose,  called  the  initial  stress  finite 
element  method  [13],  is  outlined  in  the  following. 

It  is  assumed  that  the  von  Hises  yield  criterion  [45,46]  and  flow  rules 
(See  Appendix  1)  are  valid.  Therefore,  the  following  hypothesis  appears  to 
be  generally  accepted.  If  8[e]p  denotes  the  increment  of  plastic  strain, 
then. 

*<■>,  -  x!fa}  <“> 

In  this  relation,  X  is  a  proportionality  constant  and  F  is  a  function  of 
stresses.  For  the  well-known  von  Mises  yield  surface,  F  is  given  in  the 
following  form: 

F  -  (i  («x-*y>2  +  i  <oy"«l>2  +  < V°x)2  +  3xxy  +  3xy*  +  3xlx]112  “  ®  <»>. 

where  a  is  the  effective  equivalent  stress. 

During  an  infinitesimal  increment  of  stress,  the  change  of  strain  can  be 
assumed  that  it  is  the  sum  of  the  elastic  and  plastic  parts.  Thus, 

6(e)  -  6(t]  +  6{e]  (23) 

•  P 
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(24) 


The  elastic  strain  increments  sra  related  to  stress  increments  by 
Ms}#  -  tDl"1*!®) 

Therefore,  the  strain  increment  can  be  written  as 

6(a)  -  [D]-16{cl  +  •  X  (25) 

From  the  above  equation,  it  is  convenient  to  solve  for  the  stress  changes  in 

terms  of  the  imposed  strain  changes  [13]. 

6(o]  -  ID]  6(a)  (26) 

ep 

where  [D]#p  is  the  elasto-plastie  matrix.  It  is  symmetric,  and  positive. 

«!.„  -  CO]  -  10)  CD]  [a  ♦fc&f  ID]  <27. 

where  A  is  the  hardening  parameter.  For  ideal  plasticity  with  no  hardening, 
A  will  be  xero.  If  hardening  is  considered,  in  general,  the  value  of  A  can 
be  obtained  from  the  slope  of  the  uniaxial  stress-plastic  strain  curve  at  the 
particular  value  of  plastio  strain. 

For  an  elasto-plastie  situation,  the  computation  procedures  during  a 
series  of  load  increment  can  be  outlined  as  follows: 

a)  Apply  the  load  increment  and  calculate  elastic  increments  of  stress 
[Ao’Jj.  and  the  eorr jsponding  strain  (As’}^. 

b)  Add  (Ao* },  to  the  existing  stresses  at  the  start  of  increment  {0  }  to 

x  o 

obtain  {0'). 

c)  Check  whether  F{0’}  <  0.  If  this  condition  is  satisfied,  the  pure 

elastic  case  is  continued. 
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d)  If  F{o'}  >  0  and  F{o  }  =  0,  this  means  that  the  element  was  yield  at  the 

—  o 

start  of  load  increment.  Find  {A<r}^  by  following  the  equation 
{AoK  -  [D]  tAs'}. 

1  Op  1 

e)  Evaluate  stress  which  has  to  be  supported  by  body  forces 

{Ao"^  -  {Ao'Jj  -  {Ao}1 

f)  Store  the  current  stresses  and  strains 

M  -  (o')  -  {Ao"}1 

Ce)  »  {s'}  +  ACs'Jj 

a)  If  F(o}  >  0.  but  F(o  }  <  0,  find  the  intermediate  stress  at  which  yield 

o 

begins  and  calculate  the  stress  increment  {Acr}.  Then  proceed  as  in  (d). 

h)  Calculate  nodal  forces,  or  the  plastic  load  vector  corresponding  to  the 
equilibrating  body  forces  for  any  element  by  following  relation 

{P}e  -  /[B JT{Ao/(0  jd (vol ) 

i)  Find  (Ao')2  and  (Acr'lj  using  original  elastic  properties  and  the  load 
system. 

j)  Find  the  current  value  of  a  and  repeat  steps  b)  to  i)  etc... 

The  cycling  will  be  terminated  if  the  nodal  forces  in  step  h)  reach 
sufficiently  small  values. 
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Crack 


Figure  A.  Finite  Element  Mesh 

■« 


Figure  7.  Post  Yield  (Ty  Stress  on  Crack  Plane,  y-0 

<T*  is  the  Post  Yield  Applied  Load 
v  0 
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The  J-Integral 


The  J-Integral  [8]  is  basically  a  generalization  of  Griffith's  theory  to 
elasto-plastic  rather  than  purely  elastic  crack  tip  behavior.  J  is  a  path- 
independent  integral  which  can  be  interpreted  as  the  energy  release  rate  when 
a  crack  starts  to  propagate.  In  the  finite  element  analysis  of  a  static 
crack  problem,  the  stress  intensity  factor  Kj  is  commonly  evaluated  by  the 
path-independent  J-integral  as  follows: 


;9U. 

twdy  -  Ti  sT  d3) 


where  w  is  the  elastic  strain  energy  density  and  T  is  a  path  of  arbitrarily 
chosen  contour  which  is  followed  by  the  integration  and  which  begins  at  any 
peima  on  the  lower  crack  surface,  encircles  the  crack  tip,  and  ends  at  any 
point  on  the  upper  crack  surface.  T^  is  a  traction  vector  acting  along  the 
outward  normal  to  the  contour. 


Figure  8.  J-Integral  Representation 
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SECTION  IV 


Movement  of  the  Singular  Elements 


One  of  the  latest  developments  to  numerically  model  a  propagating  crack 
is  the  use  of  moving  singular  elements.  We  treat  the  singular  element  mesh 
as  a  function  of  stresses.  When  the  Von  Mises  stress  around  the  crack  tip 
reaches  a  level  greater  than  the  yield  strength  of  the  material,  the  crack 
will  be  propagating  at  a  distance  of 


(29) 


(30) 


The  schematic  repreaentation  of  the  movement  of  the  singular  elements  is 
in  the  following  figure. 
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O  Original  Node 


#  New  Node  (After  Movement) 


Figure  9 .  Schmatic  Representation  of  the  Movement  of 
Crack-Tip  Elements 


Figure  10.  Finite  Element  Mesh  with  Quadr  ilet  erel  Quarter 
Point  Element  for  a  Center  Crack  Plate 


SECTION  V 


lUmlti  of  Nnmerical  Email 

A  rectangular  plate,  L  ■  3  in.  and  W  »  2  in.,  with  a  eantar  crack.  Fig. 
3a,  mas  considered.  The  thickness  of  the  plats  was  chosen  to  be  1/2  inch  and 
plana  stress  conditions  ware  asstuaed  in  this  ease.  The  aaterial  properties 
ware  chosen  to  be  the  following: 

E  »  10  z  106  psi 

V-  0.33 

a  •  60,000  psi 

yp 

One  to  the  symmetry  of  the  problem,  only  one-quarter  of  the  plate  needs 
to  be  analyzed.  A  total  of  163  nodes  and  44  elements  were  need  in  the 
analysis  of  growth  of  the  size  of  the  plastic  zone  with  increasing  external 
loads.  The  plastic  zones  were  plotted  by  connecting  the  centroidal  points  of 
the  yielded  elements  around  the  crack  tip. 

e 

For  the  analysis  of  the  propagation  of  the  crack,  moving  singnlar 
elements  were  need.  The  criterion  of  moving  condition  was  that  when  the 
stress  of  the  crack  tip  element  was  greater  than  the  von  Nises  effective 
stress,  the  crack  will  propagate.  The  normalized  stress  intensity  factors 
are  determined  from  averaging  the  three  different  J-integral  paths.  The 
variation  among  these  paths  is  less  than  3%. 


SECTION  VI 


~  ‘*d  Conclusions 


Although  considerable  progress  has  been  aade  in  the  field  of 
elasto-plastie  fracture  mechanics,  it  is'  the  common  consensus  of  most  of  the 
foremost  researchers  in  the  field  that  the  actual  site  and  shape  of  the 
plastic  zone  near  the  crack  tip  are  not  accurately  established  yet.  It 
would,  therefore,  be  improper  to  claim  that  the  results  presented  in  this 
report  are  the  most  reliable  and  accurate*  nevertheless,  they  provide  a  basic 
knowledge  that  can  be  used  to  better  understand  more  complex  problems  in  the 
near  future.  The  analysis  and  procedure  presented  in  this  report  for 

elasto-plastic  analysis  do  provide  a  better  and  more  reliable  method  in 
comparison  with  the  other  similar  procedures  available  in  the  literature  (12. 
17,  19,  20,  361.  For  physical  engineering  problems,  the  eight-noded 

quadrilateral  isoparametric  elements  and  the  eight-noded  quarter-point 
elements  yield  a  good  result  with  proper  calibration.  The  main  improvements 
in  using  these  elements  instead  of  the  constant  strain  triangular  elements 
are  the  convergence  and  the  stability  of  the  solution  warranted. 

It  has  been  demonstrated  that  the  technique  can  be  applied  where  crack 
propagation  is  very  small  and  may  not  otherwise  be  measured  by  conventional 
experimental  techniques. 

The  following  conclusions  are  based  on  the  elastic-plastic  crack  growth 
analysis. 

1)  A  method  was  developed  to  calculate  the  size  and  shape  of  the  plastic 

zone  around  the  crack  tip.  The  problem  was  solved  by  using  the  initial 

stress  finite  element  method.  Eight-noded  quadrilateral  isoparametric 
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•  1  wants  war*  usad  in  tha  analysis.  It  was  provsd  that  this  nathod  is 
quit a  affieiant  and  provided  raasonabla,  accnrata  rasnlts. 

2)  Tha  eight-noded  qoartar-point  element  yialdad  a  good  rasult  with  propar 
oalibration.  This  singular  element,  whan  coupled  with  the  moving 
a  1  asiant  technique,  was  praf  arable  for  studying  the  fast  fracture 
problem. 

3)  It  was  found  that  tha  procedure  of  tha  moving  singular  element  was  much 
easier  to  use  than  tha  so-called  crack  tip  node  release  method. 
Furthermore,  tha  moving  singular  alamant  provided  realistic  crack  growth 
values. 

4)  Tha  moving  finite  alamant  technique  can  be  modified  to  solve  the  crack 
branching  problem. 

Tha  above  advancements  in  tha  understanding  of  elastic-plastic  crack  growth 
are  specially  suited  for  aiding  tha  future  study  of  tha  fast  or  dynamic 
fracture  problems.  In  addition,  the  crack  growth  criteria  investigation  and 
tha  calculation  of  the  plastic  zona  size  provide  significant  progress  towards 
life  prediction  of  machine  elements,  such  as  turbine  blades  and  turbine 
di sks . 

The  following  additional  research  work  is  recommended  to  further  the 
machine  element  life  prediction  capability. 

1)  Future  work  needs  to  include  impact  loading  and  cyclic  loading 
conditions. 

2)  Use  present  approach  to  analyze  other  test  specimen  geometries  to 
determine  the  dependence/ independence  of  results  on  speciment  geometry. 
Also,  investigate  the  repeatability. 
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3)  Um  present  approach  to  analyze  the  crack  growth  of  the  composite 
■aterial  structures. 

4)  Environmental  effects  should  be  further  researched  by  conducting  fatigue 
and  creep  crack  growth  tests  in  different  environmental  conditions,  such 
as  high  sulfur  content,  vacuum,  etc. 

In  summary,  the  technique  developed  here  is  worth  exploring  further  to 
provide  a  better  understanding  of  crack  growth  rate  behavior  in  machine 
elements  under  normal  operating  conditions. 

SECTION  VII. 


Appendix  I 


,  1.  Von  Ml see  Criteron 

A  criterion  that  involves  a  smooth  function  was  proposed  by  Von  Hises. 
In  its  most  widely  used  form,  the  von  Mises  criterion,  in  terms  of  principal 
I  stresses,  prediets  that  yielding  occurs  when 

Ss 

(oj  -  a 2>2  +  (<t2  -  Oj)2  +  (»3  -  Oj)  ■  constant  (1A) 

I  In  a  more  general  form, 

(<»  -  a)2  +  (o„  -  a)2  +  (o„  -  a)2  +  6(t  2  +  x  2  +  x  2)  “  constant  (2A) 

x  y  y  z  z  x  xy  yz  zx 

To  determine  the  constant,  the  following  procedure  is  used: 

a)  For  uniaxial  tension,  yielding  occurs  when  o.  *  o _ ,  o,  ■  a.  *  0. 

1  yp  2  3 

Therefore  equation  (1A)  becomes 

2  2 
*"»  2®,  ■  constant  •  2o _ 

r  1  ” 

» 

li 
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•  0  and  using 


For  pur*  shear,  yielding  occurs  when  e,  ■  -  a,  •  k,  «2 
equation  (1A): 

2  A  2  A  .  2  .  2  4  t  -.2 

Oj  +  Oj  +  40^  ■  0Oj  ■  constant  ■  6k 

Thus  the  von  Miscs  criterion  nay  be  written  as  : 

(Oj  -  o2)2  ♦  <o2  -  ®3 > 2  +  <*3  -  ffj)2  "  2oyp2  -  <k2 


According  to  this  criterion  the  tensile  and  shear  yield  stresses  are  related 
as  o  ■  J3k. 

It  is  convenient  to  consider  this  criterion  as  a  function  of  an 
effective  stress,  denoted  as  a",  where  o  is  a  function  of  the  applied 
stresses.  Whenever  its  magnitude  reaches  the  yield  strength  in  uniaxial 

tension,  then  that  applied  stress  should  cause  yielding  to  occur.  Thus 

—Ir/  .2  ^  .  .2  .2.1/2  .... 

a  •  “  a2^  *  +  ®3  “  °i )  1  (4A) 

I*3 

V1* 


<ryi 


!  Ap 


•  Jj  fur.  Shear 


fcr„  •  <r, 


T  *<- »r»:  ■<£, 


P i|«r.  IA.  Von  HUos  Ttold  Locum  Oj*0 


An  interpretation  of  the  Hises  criterion  is  that  yielding  occurs  when 
the  elastic  energy  causing  distortion  reaches  a  critical  value.  This  strain 
energy  is  found  in  a  general  way  by  subtracting  the  dilatational  strain 
energy  from  the  total  elastic  strain  energy.  The  total  strain  energy  per 
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unit  volume  is: 

D  »  1/2  (o  8  +08  +oe  +t  v  +  t  y  +  t  y  ) 

v  z  x  y  y  z  z  xy  xy  yz  yz  zz  zz 


or  for  the  esse  of  principsl  stresses: 


ny  -  1/2  (o^j  +  o2e2  +  o3s3) 


(SA) 


(6A) 


To  express  the  shove  eqastion  ss  s  function  of  stresses,  the  generslized  form 
of  Hooke's  lev,  gives 

0  \) 

(7A) 


1  2  2  2  V 

Uy  «  2E  +  ®2  +  )  — g  ^°\a2  *  a2a3  + 


Since  only  normsl  stresses  esnse  s  volume  change,  the  dilststion  is: 
A  -  8i  +  e2  +  83“  v g—  (ox  +  o2  +  o3)  »  -|  (l-2v)  om 


(8A) 


Let  A  *  3e 

m 

l-2v 

Therefore  8  *  ~ s —  O 

m  E  m 

From  the  shove  discussion,  it  esn  he  seen  that  the  work  due  to  dilststion  is 

0.  -  3(1/2  oe). 
a  m  m 


then  U.  -  t3(l-2V  )a  ]/2E 
d  a 


l-2v 


d  6E  '°1  +  a2  +  V 


(9A) 


By  subtracting  eq.  (9A)  from  (7A)  to  give  the  shear  strain  energy,  U#,  the 


result  after  rearrangement  is: 

Us  ■  T2B  I(<T1  “  ff2) 2  +  (c2  ■  °3)2  +  (ff3  "  ffl)2] 


(10A) 


Nov,  the  shear  strain  energy  induced  during  uniaxial  tension  where  o2  *  o, 
0  is: 
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u  - 

s 


The  critical  value,  U  »  that  may  be  developed  to  cause  yielding  will  result 

SC 

when  0,  *  a  .  Setting  eq.  (10A)  equal  to  this  critical  value  leads  to: 

1  yp 


t(Oi  -  a2)2  +  <«2  -  o3)2  +  (o3  -  Bl)2]  = 


(11  A) 


3 .  Flow  Rules  or  Plastic  Stress-Strain  Relationships 


Consider  plastic  flow  under  uniaxial  tension  as  indicated  in  following 


figure: 


«l£i 


01-0. 

ae^o 


<r}=o,  a£j%  0 


Figure  2A.  Stresses  and  Incremental  Plastic  Strains 
for  Uniaxial  Tension 


Now  the  deviatoric  stress  in  the  1  direction  is  0^  ■  0^ 

particular  instant  represented  in  the  above  figure. 


-  0  and  at  the 


0'j  ■  2/3  0j  and  o^' 


0-1/3  Oj  -  -  1/3  0X 


Therefore, 


-  2o2'  -  -  203' 
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volume  constancy,  the  sume  of  the  plastic  strain  increments  must  be 
>ref  ore, 

+  de2  +  de3  =  0 

:re  is  symmetry  in  this  instance,  de^  =  d e „ ,  therefore,  de^  =  -  2de2 
This  leads  to: 


constant 


dX 


(12a) 


ication  is  that  the  ratio  of  the  current  incremental  plastic  strain 
s  to  the  current  deviatoric  stresses  is  a  constant.  This  expresses 

t 

dtl-Iteuss  flow  rules,  where  the  elastic  strain  increments  have  been 

convenience,  the  flow  rules  may  be  expressed  in  other  forms 

^  -a, 

1  -  °2 

*3  dA[0l  -  2  <02  +<’3>1  (13a) 

■  ^  1°1  -  7  <«2  +  ”3)1 


is  the  incremental  effective  strain  which  is  defined  as 

n  2  9  ,1/2 

“  -  d£2>  +  (d£2  -  de3>  +  (de3  -  dei>  ] 


(14a) 
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